Astronomy & Astrophysics manuscript no. clustmass8 


© ESO 2010 


September 6, 2010 





Calibration of radii and masses of open clusters with a simulation. 

A. Ernst 1 , A. Just 1 , P. Berczik 1 ' 2 ' 3 , and M. I. Petrov 4 

1 Astronomisches Rechen-Institut, Zentrum fiir Astronomie der Universitat Heidelberg, Monchhofstrasse 12-14, D- 
69120 Heidelberg, Germany 

2 National Astronomical Observatories of China, Chinese Academy of Sciences, Datun Lu 20A, Chaoyang District, 
Beijing 100012, China 

3 Main Astronomical Observatory, National Academy of Sciences of Ukraine, Akademika Zabolotnoho 27, 03680 Kyiv, 
Ukraine 

4 Institut fiir Astronomie der Universitat Wien, TiirkenschanzstraBe 17, A- 1180 Wien, Austria 
Received ... Accepted ... 

ABSTRACT 

Context. A recent new approach to apply a simple dynamical mass estimate of tidally limited star clusters is based on 
the identification of the tidal radius in a King profile with the dynamical Jacobi radius. The application to an unbiased 
open cluster catalogue yields significantly higher cluster masses compared to the classical methods. 
Aims. We quantify the bias in the mass determination as function of projection direction and cluster age by analysing 
a simulated star cluster. 

Methods. We use direct Af-body simulations of a star cluster including stellar evolution in an analytic Milky Way 
potential and apply a best fit to the projected number density of cluster stars. 

Results. We obtain significantly overestimated star cluster masses which depend strongly on the viewing direction. The 
overestimation is typically in the range of 10-50 percent and reaches a factor of 3.5 for young clusters. Mass segregation 
reduces the derived limiting radii systematically. 

Key words. Galaxy: open clusters and associations: general - methods: Af-body simulations - stellar dynamics 



1. Introduction 

In a series of papers (Piskunov et al. 2007, 2008a, b) a new 
approach to determine the masses of open star clusters 
(OCs) was developed and used to determine the initial and 
present day mass function of OCs in the solar neighbour- 
hood. The new method is based on the determination of 
the tidal radius r t from the cumulative number of cluster 
members as function of projected distance to the cluster 
center. For each cluster the tidal radius r t is determined 
from projected number density profiles by fitting a King 
1962 profile (King 1962). The identification of the King cut- 
off radius r t with the "Jacobi" radius r j (i.e. the dynamical 
tidal radius, which is the distance from the cluster center to 
the Lagrange points L\ and L<i) yields then the OC mass 
from the standard formula (Equation [9] below solved for 
M c \). The application of this dynamical mass estimate of 
tidally limited clusters to an unbiased OC catalogue yields 
an independent mass determination compared to the classi- 
cal methods. A detailed comparison with other methods of 
cluster mass determinations is also given. In a second step 
the method is extended to all OCs of an unbiased cluster 
catalogue by establishing a transformation of the observed 
semi-major axis and central surface density to r t . These re- 
sults were then used to derive the cluster present day mass 
function (CPDMF) and the initial mass function of OCs 
(CIMF) in the extended solar neighbourhood. Adopting a 
constant cluster formation rate over the last 10 Gyr yields 
a surface density of 18M0pc -2 of stars born in OCs. This 
corresponds to a fraction of 37% of disc stars which were 
born in OCs (Roser et al. 2010). This is large compared to 



the classical values of the order of 10% or less (e.g. Wielen 
1971, Miller & Scalo 1978). 

Some crucial assumptions enter the dynamical mass de- 
termination based on fitting a King profile: a) The OC fills 
its Roche lobe in the tidal field of the Milky Way. For 
compact (e.g. Roche-lobe underfilling) clusters rj and as 
a consequence the mass can be underestimated by a large 
amount, b) The effect of mass segregation can be neglected, 
i.e. star counts of the upper main sequence, which domi- 
nate the observed cluster members, are representative for 
the mass distribution, c) The elliptic shape of the clusters 
and the contamination through tidal tail stars do not re- 
sult in a systematic bias. Shape parameters were measured 
by Kharchenko et al. (2009) and the distribution of tidal 
tail stars were investigated in detail (e.g. Just et al. 2009). 
d) The tidal radius r t determined by fitting the cumula- 
tive projected mass profile represents the Jacobi radius rj 
to derive the cluster mass. Since the cluster mass depends 
on the third power of rj, the method is very sensitive to 
systematic errors in the derivation of rj. 

In the present paper we quantify the possible bias in- 
troduced by the identification of the tidal radius from a 
King profile fitting r t with the Jacobi radius rj used for 
the mass determination by Piskunov et al. (2007). We ap- 
ply the King profile fitting procedure to a direct A/"-body 
simulation of a dissolving star cluster at different evolution- 
ary states. We have simulated a star cluster on a circular 
orbit at Rc = 8.5 kpc which evolved in the tidal field of the 
Milky Way including stellar evolution. We took snapshots 
of the evolved model with all stellar masses and positions 
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Table 1. The list of galaxy component parameters. The 
first column gives the component, the second the mass, 
and the third and fourth the Plummer-Kuzmin parameters 
(equation pi) . 



Component 


M [M ] 


a [kpc] 


b [kpc] 


Bulge 


1.4 x 10 10 


0.0 


0.3 


Disk 


9.0 x 10 10 


3.3 


0.3 


Halo 


7.0 x 10 11 


0.0 


25.0 



at four different times and projected the snapshots from 
the perspective of an observer on Earth (at Ro = 8 kpc) 
onto the sky, at different positions along its orbit. After all, 
we determined the model's limiting radius r t by fitting the 
projected cumulative mass profile with Equation (|6| and 
compared r t to the actual Jacobi radius rj. 

The paper is organized as follows: In Section 2 we dis- 
cuss the method of TV-body simulations in the external po- 
tential of the Milky Way, Section 3 contains the theory of 
the cluster geometry in a tidal field. In Section 4, we show 
a simple iterative method to determine the Jacobi radius of 
an Af-body model of a star cluster in a tidal field. In Section 
5 the projection and fitting methods are described. Finally, 
Section 6 contains the results and section 7 the conclusions. 
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Fig. 1. Top: Rotation curve (at z = 0) of the 3- 
component Plummer-Kuzmin model of the Milky Way. 
Bottom: Epicyclic and vertical frequency parameters j3 = 
k/CI and 5 = u/fl (at z = 0). 



2. Numerical simulation 

We analyse in detail a numerical simulation of a star clus- 
ter with initial mass Mq = 1O 4 M0 on a circular orbit in an 
analytic Milky Way potential. It is the fiducial cluster sim- 
ulation (run 10) discussed in Just et al. (2009). We have 
chosen this cluster, because it is a typical representative 
for the high-mass end of the observed OCs. Since the total 
mass of the cluster system is dominated by the high-mass 
end, the correction of biases in the mass determination are 
most important in that parameter regime. The cluster is 
set up as a Wo = 6 King model with a half-mass radius of 
8 pc. The extension of the cluster initially exceeds the Roche 
lobe initially leading to an enhanced mass loss in the first 
0.5 Gyr. We used a Salpeter IMF and included mass loss 
by stellar evolution. The total lifetime of the cluster at a 
circular orbit with Rq = 8.5 kpc is 6.3 Gyr. For details of 
the evolution see Just et al. (2009). 

For the high-resolution simulation of a dissolving star 
cluster with N = 40404 particles in the tidal field of the 
Milky Way the direct TV-body code ^grapeQ (Harfst et al. 
2007) has been used in combination wit the micro-GRAPE6 
special-purpose hardware at the Astronomisches Rechen- 
Institut (ARI) in Heidelberg] (/>grape is an acronym for 
Parallel Hermite Integration with grape. The code is writ- 
ten in ANSI-C and uses a fourth-order Hermite scheme 
(Makino & Aarseth 1992) for the orbit integration. It is par- 
allelized and uses the MPI library for communication be- 
tween the processors. The force computations are executed 



1 The present version of the code is pub- 
licly available from one of the authors FTP site: 
f tp : / /ftp . ari . uni-heidelberg . de/ staff /berczik/ 
phi-GRAPE-cluster/code-paperA 



on the fast special-purpose hardware grape. The special- 
purpose micro-GRAPE6 hardware cards are especially de- 
signed to calculate gravitational forces in iV-body simu- 
lations very fast using parallelization with pipelining (see 
Harfst et al. 2007 and references therein). 

The code 0GRAPE does not use regularization as 
the codes nbody4 or nbody6++ (Aarseth 1999, 2003; 
Spurzem 1999) but a standard Plummer type Af-body grav- 
itational softening. The softening length in the model used 
for the current work was e = 10 -3 pc. We tested with differ- 
ent softening lengths e = 10 -3 , 10 -4 and 10 -5 pc that there 
are no significant differences regarding shape evolution and 
star cluster mass loss. 

For the simulation of a star cluster in the tidal field 
of the Galaxy the A"-body problem is solved in an ana- 
lytic background potential. We use an axi- symmetric 3- 
component model, where bulge, disc and halo are described 
by Plummer-Kuzmin models (Miyamoto & Nagai 1975) 
with the potential 



$(R,z) 



GM 



R 2 + (a + Vb 2 + z 2 ) 2 



(1) 



GRACE: http : //www . ari . uni-heidelberg . de/grace 



The parameters a, b and M of the Milky Way model are 
given in Table [T] for the three components. 

The top panel of Figure [I] shows the rotation curve of 
the 3-component model of the Milky Way. The parameters 
of the 3-component model are chosen such that the rotation 
curve matches that of the Milky Way (Dauphole & Colin 
1995). At the solar radius Rq = 8.0 kpc, which was assumed 
in this study, the value of the circular velocity is Vq = 234.2 
km/s. The values of Oort's constants A and B are consistent 
with the observed values (A, B) = (14.5 ± 0.8, -13.0 ± 1.1) 
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km/s/kpc derived by Piskunov et al. (2006). More gener- 
ally, the dimensionless epicyclic and vertical frequency pa- 
rameters are given by 



l.l 



(3 2 = k 2 IVl 2 = 2 

s 2 = u 2 /n 2 



d\nn 
d\nR 



and 



n 2 



(2) 
(3) 



where ft, v and Q are the epicyclic, vertical and circular fre- 
quencies of a near-circular orbits and p is the local galactic 
density (see Oort 1965 for the derivation of 5 2 ). The bot- 
tom panel of Figure [I] shows the course of the epicyclic and 
vertical frequency parameters f3 and 5. 

At the radius Rc = 8.5 kpc of the circular orbit consid- 
ered in this study we obtain (/3c, Sc) = (1-37, 2.86) and the 
circular velocity Vq = 233.3 km/s. The orbital time scale 
at R c = 8.5 kpc is T orb « 224 Myr. 



3. Cluster geometry 

According to King (1962), the projected density profile T>(r) 
of a star cluster can be approximated by 

E(r) = /c{x- 1/2 -C- 1/2 } 2 for r < r t (4) 

with normalisation constant k and 

X(r, r c ) = l + (r/r c ) 2 and C(r c , r t ) = 1 + (n/r c ) 2 , (5) 

where r c is the core radius and r t is the (tidal) cutoff radius 
where the projected density of the model drops to zero. 
Integration yields the cumulative form of the King 1962 
profile, 



M p (r) = 2n [ E(r')r'dr' for 
Jo 



r 2 k 



ln(X) 



r <r t 



1 X 



C l/2 



c 



(6) 



where M p (r) is the projected cumulative mass of the model, 
i.e. the mass in projection on the sky within a circle of ra- 
dius r. For r > r t we force the integrated profile (Equation 
|6| to the finite value 



M p (r t ) 



<- 2 k 



ln(C7) 



C 1 / 2 C 



(7) 



We will not use the identification of M p (r t ) with the clus- 
ter mass, since in practice the equations are applied to star 
counts and the effective mass-of-light ratio enters the nor- 
malisation constant k. 

Already King (1961) remarks that the tidal forces from 
the galaxy distort only the outer regions of a star cluster. 
We quantify these deviations from spherical symmetry due 
to the tidal field in this Section. We employ a coordinate 
system (x,y,z) of "principal axes of the star cluster". Its 
origin is the cluster centre. The x-axis points away from 
the galactic centre, the y-axis points in the direction of 
the galactic rotation and the z-axis is directed towards the 
galactic north pole. Figure [2] shows a "principal axis plane 
cut" through the equipotential surfaces of the effective po- 
tential around a star cluster on a circular orbit in the tidal 
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Fig. 3. Principal axis ratios of the equipotential surfaces 
around the cluster as a function of the parameter 7 = x/rj. 
At the centre of the cluster we have 7 = while 7 = 1 
(dashed line) corresponds to the critical equipotential sur- 
face through x = rj. We assumed (/?, S) = (1.37,2.86) and 
a Kepler potential for the cluster. 



field of the Milky Way. To second order, the effective po- 
tential is given by 



eff 



$eff,0 



GM C] 



\/x 2 +y 2 

2 - 4)n 2 x 2 - 



s 2 n 2 z 2 



(8) 



For the cluster we assumed a Kepler potential, which is 
a very good approximation in the outer parts (Just et al. 
2009). The unit in Figure [2] is the Jacobi radius Tj. The 
Jacobi radius is defined as the distance from the cluster 
centre to the Lagrange points L\ and L2. It is given by 



rj 



GM C 



(4 - /3 2 )n 2 



1/3 



(9) 



(see King 1962). The value of the effective potential on the 
critical equipotential surface which connects L\ and L2 can 
be easily calculated from (|8| and (|9|. It is given by 



^eff,crit = $eff,0 



3GM C] 



2 rj 



(10) 



Assuming that x max = rj, the radius in ^/-direction of the 
last closed (critical) equipotential surface follows from (JsJ) 
and JTob, 



(11) 



For the radius in z-direction we have to solve a cubic equa- 
tion 



K 



with the constants 



- 1-22 _ 

) ^ ^max — 



3 if 
2M 



K = GM ch L = SQ, M = rj. 
The only real solution is given by 



(12) 



(13) 
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{K [K + 2L (LM 3 + VKM 3 + L 2 M 6 )} } 1/3 - K 



2/3 



L (LM 3 + Vif M 3 + L 2 M 6 ) 

For (/3, J) = (Pc,6c) = (1.37,2.86) we obtain 
^max « 0.503 r j 



1/3 



(14) 



(15) 



(e.g. Wielen 1974). More generally, the ratios of the princi- 
pal axes ?/ max : x max and z max : x max of all closed equipo- 
tential surfaces are shown in Figure [3] as a function of the 
parameter 7 = x/rj. At the centre of the cluster we have 
7 = while 7 = 1 corresponds to the critical equipotential 
surface through x = rj. 

Figure [4] shows the logarithmically colour-coded surface 
density E of projections of the simulated A'-body model of 
the star cluster onto its principal axis planes at time T4 = 
1.31 Gyr. The surface density E has been calculated from 
the A'-body snapshot file with the method of Casertano & 
Hut (1985) (their Equation II.6 with j=20). The contours 
correspond to A log E « 2 dex. The extent of the last closed 
(critical) equipotential surface is marked with dashed lines. 
The contours of constant surface density roughly follow the 
equipotential surfaces from Figure [2] 

4. Jacobi radius of the model 

We determined the Jacobi radius of our simulated A'-body 
model iteratively from Equation [9j The gravitational con- 
stant G and the quantities j3 and Vt are known but M c i, the 
cluster mass within r = Tj is unknown since r j is unknown. 
Our iteration is given by 



'n+l 



GM enc (r n ) 
_(4-/3 2 )^ 2 



1/3 



(16) 



where M enc (r) is the enclosed mass of cluster stars within 
radius r around the cluster centre. Starting with ro = 00 
the iteration theoretically converges towards = rj. In 
practice, a few iterations are sufficient to determine rj ac- 
curately. It is interesting to note that this simple method 



does not rely on numerical fits of the effective potential to 
the envelope of the spatial distribution of Jacobi energies 
of the cluster stars. It can easily be applied to an A^-body 
snapshot file which contains masses and positions of the 
individual particles at a certain time. 



5. Projection and fitting 

We calculate the projection at the sky in Galactic coor- 
dinates (Z,6) and ignore here for simplicity the perspective 
effects. Strictly speaking the projection of the cluster model 
is consistent only in the midplane 6 = 0° corresponding to 
the circular orbit in the Galactic plane. The rotation angle 
a is related to the galactic longitude I by the law of sines, 



R 

sin a = — sin I — — Rc > Ro 
Rc 



(17) 



with ^0 = 8 kpc for the solar circle and Rc = 8.5 kpc for 
the cluster orbit. At the Galactic coordinates (l,b) = (0,0) 
and (l,b)= (180°, 0) the "system (x,y,z) of principal axes of 
the cluster" with origin at the cluster centre is not rotated. 
Note that at I = 90° and I = 270° we have the maximum 
rotation angle a = —70°. 25 and a = 70°. 25, respectively. 

A sketch of the projection of the simulated A^-body 
model of the star cluster onto the plane of the sky is illus- 
trated in Figure |5j Four orbital positions of the dissolving 
star cluster with its tidal tails are marked in the sketch. 

In order to demonstrate the effect of different projec- 
tions also in galactic latitude, we add the cases of a certain 
height of the cluster orbit above and below the plane of 
the solar circle (dotted lines in Figure [5]). The effect of the 
corresponding vertical oscillation of the cluster orbit on the 
intrinsic structure of the cluster and of a variation of S in 
Equation [3] are neglected in this study. 

For all positions of the cluster on its orbit, we rotate the 
cluster around the z axis and then around the y' axis by 
the ordered pair of angles (a, -b) in order to simulate the 
perspective of the cluster for an observer on earth. After 
the projection, we determine the polar symmetric profile 
of the projected cumulative mass M p (r) from our A^-body 
data file by summations over radius and polar angle and 
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Fig. 4. Surface density of projections onto the principal axis planes of the cluster at T4 = 1.31 Gyr. The dashed lines 
show the theoretical values of x max /rj, y max /rj and £ max /rj. The contours correspond to AlogS ~ 2 dex. 




Fig. 5. Sketch of the galactic coordinate system in which 
the projection is done. The galactocentric radius of the sun 
is assumed to be Ro = 8.0 kpc while the cluster orbits at 
Rc = 8.5 kpc on a circular orbit. The abbreviations denote 
the Galactic North (NGP) and South (SGP) Pole and the 
Galactic Centre (GC). 



apply a fit with Equation [6] For the fitting, we used the 
mpfit package in idl (Markwardt 2009; More 1978 for the 
Levenberg-Marquardt algorithm) . 



6. Results 

We first discuss projections according to Equation [TT] in 
detail. Figure [6] shows examples of fits (upper panels) with 
the corresponding projections (lower panels). The resulting 
parameter ratio r t /rj is given in the upper panels and the 
(l,b) coordinates in the corresponding lower panels. In the 
upper panels, the solid (black) line represents the data and 
dotted (blue) line the fit. The dashed (red) lines mark r c 
(left dashed line) and r t (right dashed line) from the fit 
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Fig. 8. Histogram of the fraction on the sky per bin in 
Mt/M c i corresponding to the T4 = 1.31 Gyr parameter 
surface from Figure [7| M t and M c \ are the "tidal masses" 
calculated with Equation [9] from r t and rj, respectively. 




Fig. 9. Mean mass (M t /M c \) as a function of the galactic 
latitude b. M t and M c \ are the "tidal masses" calculated 
with Equation |9] from r t and rj, respectively. 



with Equation [6] In the lower panels, the dashed (red) line 
marks Tj. 
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Fig. 6. Examples of fits (upper panels) with the corresponding projections (lower panels) at time T4 = 1.31 Gyr . The 
resulting parameter ratio r t /rj is given in the upper panels and the (l,b) coordinates in the corresponding lower panels. 
In the upper panels, the solid (black) line represents the data and dotted (blue) line the fit. The dashed (red) lines mark 
r c (left dashed line) and r t (right dashed line) from the fit with Equation [6j In the lower panels, the dashed (red) line 
marks rj. 
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Fig. 7. Parameter surfaces of r t /rj as a function of Galactic coordinates for a fit of King (1962) models to projections on 
the sky of a simulated model at different positions on its theoretical orbit. We used a squeezed Hammer- Aitoff projection. 
The color denotes the value of r t /rj on a linear scale. The plots in the top row correspond to T\ = 0.62 Gyr (top left) 
and T2 = 0.84 Gyr (top right). The plots in the bottom row correspond to T3 = 1.06 Gyr (bottom left) and T4 = 1.31 
Gyr (bottom right). 



We derived the ratio r t /rj for all projection directions 
in (Z, b) at four different evolution times of the cluster. The 
projections are done in steps of 4 degrees from b — —90° 
to b = +90° and Z = to I = 360°. The full TV-body 
snapshot file with N = 40404 particles has been used for the 
projection and fitting procedure as described in Section 5. 
Figure [7] shows parameter surfaces of r t /rj as a function of 
Galactic coordinates of the cluster centre at times T\ =0.62 
Gyr (top left), T 2 = 0.84 Gyr (top right), T 3 = 1.06 Gyr 
(bottom left) and T4 = 1.31 Gyr (bottom right). A squeezed 
Hammer- Aitoff projection of Galactic coordinates has been 
used (see Appendix A). The total number of fits for each 
plot in Figure [7] is Nf[ ts = 4186. Note that in the plo t for 
Ti, a peak around (Z,6) w (270°, 0) (see also Figure A.l 
in Appendix A) is not resolved (white colored area). Here 
r t /rj reaches a factor of 1.52. 



It can be seen that the cluster masses are typically over- 
estimated. The strongest bias to high masses at time Ti 
(including the peak) shows that the unbound stars stay for 
a long time close to the cluster and contribute to the outer 
density profile in the fitting procedure. This effect depends 
on the initial conditions. At later times T3, T4 the cluster 



mass distribution becomes stable and the bias is indepen- 
dent of the age of the cluster. 

Figure [8] shows a histogram of the fraction on the sky 
per bin in M t /M c i = (r t /rj) 3 for the parameter surface 
in Figure [7] corresponding to time T4, where the "tidal" 
masses M t and M c \ are calculated with Equation [9] from r t 
and rj, respectively. The solid line shows the distribution 
for projections in the range between —20° < b < +20°. The 
most frequent overestimation of the mass is M t /M c i = 1.18. 
The dashed line shows the distribution for projections in 
the range between —40° < b < +40°. The most frequent 
ratio is the same as for the solid line. However, the ratios 
extend up to M t /M c i = 1.7. The dotted line shows the 
distribution for all projections. The most frequent ratio is 
the same as for the other two lines, but the ratios extend 
up to Mt/Md « 2.1. 

Figure |9] shows the mean mass (M t /M c \) as a function 
of the galactic latitude 6, where M t is the mass of cluster 
stars within radius r t which has been obtained by the fitting 
procedure. Shown are the curves corresponding to times 
Ti — T4. All curves show a local minimum which is located 
roughly at 6 = 0°. For time T\ the mean overestimation of 
the mass reaches 2.5 in the direction of the Galactic poles. 
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Fig. 10. Parameter surfaces of r t /rj and Ar t /r t as a func- 
tion of Galactic coordinates for a fit of King (1962) models 
to to projections on the sky of a simulated model at differ- 
ent positions on its theoretical orbit. We used a squeezed 
Hammer- Aitoff projection. The time is T4 = 1.31 Gyr. The 
color denotes the value of r t /rj (and Ar t /r t , bottom plot) 
on a linear scale. The upper plots shows the parameter 
surface for the 400 brightest stars in the simulated cluster 
(data courtesy of J. Beuria). The middle and bottom plots 
show the result of the bootstrap analysis (for explanations 
see the text). The middle plot shows the mean value of 
r t /rj averaged over 100 small- number samples. The bot- 
tom plot shows the corresponding relative standard devia- 
tion Art/rt. 



Usually only a small fraction of cluster stars are iden- 
tified as members leading to an increased statistical un- 
certainty in the fitting procedure. The uppermost plot in 
Figure [10] shows the parameter surface derived for the 400 
most massive (i.e. the brightest) stars in the simulation at 



time T4 (data courtesy of J. Beuria). The reason for the low 
values of r t /rj is mass segregation. The stronger asymme- 
tries in ±b compared to the plots in Figure [7] are due to a 
slightly asymmetric distribution of the 400 mass-segregated 
stars in position space. 

In order to measure the statistical scatter in r t we ap- 
plied a bootstrap analysis. We divided the TV-body snap- 
shot file for time T 4 with N = 40404 particles into 100 
small-number samples of Af samp i e = 400 particles each (the 
remaining particles were dropped out of the analysis). For 
each of these small-number samples we applied the proce- 
dure described above. The resulting total number of pro- 
jections and fits was therefore iVfits = 418600. The mid- 
dle and bottom plots in Figure 10 show the result of the 
bootstrap analysis. The middle plot of Figure [10| shows the 
mean value of r t /rj averaged over the 100 samples. The 
bottom plot shows the relative standard deviation Ar t /r t 
which resulted from the averaging over the 100 samples. 
The uncertainty on a single determination of r t lies in the 
range between 10 and 20 percent. The highest uncertainty 
is expected for redeterminations in the vicinity of the peak 
around (Z,6) « (270°, 0). 



A comparison with the uppermost plot in Figure 10 
shows that the derived limiting radii for the sample of the 
most massive stars are systematically lower due to mass 
segregation. The differences are typically on a 1 — 2a level. 

7. Conclusions 

The result of the analysis in this paper is the confirmation 
that the star cluster masses are typically overestimated if 
the method by Piskunov et al. (2007) is applied on a com- 
plete sample. Moreover, we quantified the methodological 
error in our analysis. 

Figures [7] and [10] show that at certain Galactic coordi- 
nates the King (1962) profile fits are particularly biased. A 
high bi as is predicted for (Z, b) « (270°, 0) (see Figures [7| 
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and A.l). The corresponding rotation angle of the clus- 
ter is a = 70°, where the projection is parallel to the inner 
end of the tidal arms (see lower left plots in Figure [6]). For 
(Z, b) = (90°, 0) there is no corresponding peak in the pa- 
rameter surface (best visible for T\) due to the asymmetry 
between the leading and trailing tidal tails (see the sketch 
in Figure [5]). Al so, for Galactic latitudes beyond b ~ ±40° 
(cf. Figure~ [A?l| the bias becomes large but only few OCs 
are located in that regime. 

For the parameter surfaces in Figure [7] the masses are 
biased within the ranges [1.3 M ch 3.5 M d ] (at Ti), [1.1 M ch 
2.3 M c i] (at T 2 ), [1.0 M ch 2.0 M cl ] (at T 3 ) and [1.0 M ch 
2.1 M c i] (at T4) depending on the evolutionary state of the 
star cluster in the tidal field of the Galaxy. The bias de- 
pends strongly on the projection angles, which transform 
differently to Galactic coordinates for different orbital radii 
of the OC. 

Furthermore, a bootstrap analysis showed that the rel- 
ative error on a single determination of the limiting radius 
r t lies in the range between 10 and 20 percent (at time 
T4) corresponding to an uncertainty in the mass of ~ 50 
percent for samples of N samp \ e = 400 particles (which are 
typical for rich OCs). 

Mass segregation of the brightest stars in a cluster can 
alter the r t /rj factor significantly which is important for 
the data analysis of observations. The mass segregation re- 
sults in a concentrated core which leads to an underesti- 
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mation of the tidal radius. For a younger cluster age one 
would expect lesser mass segregation. 

For a quantitative correction of the bias in the cluster 
mass determination by identifying r t with r j an extensive 
parameter study of cluster parameters is necessary. The 
influence of mass-segregation on selection effects concerning 
the brightness limit of the observations should be included, 
because stellar evolution is taken into account. Especially 
for young clusters the bias factor can be sensitive to the 
initial conditions. The final goal is to find an agreement 
of OC mass determinations by the different methods. This 
would allow an interesting insight in the OC properties like 
the IMF, mass-to-light ratio and mass segregation. 
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Fig. A.l. Squeezed Hammer- Aitoff projection. The 
Galactic latitude runs from b = —90° to b = +90° and the 
longitude from I = to I = 360°. 

Appendix A: Squeezed Hammer-Aitoff projection 

The squeezed Hammer-Aitoff projection is given by 



x = 2f 



cos(6) sin(//2) 



y/l + cos(6) cos(//2) ' 
2 sin(6) 

7 y/l + C0S(6) COs(//2) ' 



(A.l) 
(A.2) 



It is the standard Hammer-Aitoff equal-area projection 
where we introduced a free squeezing factor /. The squeez- 
ing leaves the area element dA = dxdy = y/detg dbdl = 
cos(b)dbdl invariant, where g is the fi rst f undamental form 
calculated from equations (A.l) and (A.2). The ratio of di- 



ameters of the elliptic projection area is given by d x /d y — 
f 2 . In the standard Hammer-Aitoff projection we have 
/ = y/2. For a projection onto a circular area one can set 
/ = 1. The inverse projection is given by 



x 

V 

arcsm [ {^zy 



fy 

2 



b 

I = 2 arctan 



1 zx 
2fz 2 -l ) ' 



(A.3) 
(A.4) 
(A.5) 



where we have introduced the auxiliary variable z. 

We used / = 1.2 for Figures [7] and [To] The resulting 
coordinate system (which is hidden in Figures |7|and[l0) can 
be seen in Figure A.l We have modified idl routines by 
W. B. Landsman to incorporate the free squeezing factor. 
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